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Abstract We propose an analytical model to predict the adhesives pene- 
tration into hard wood. Penetration of hard wood is dominated by the vessel 
network which prohibits porous medium approximations. Our model considers 
two scales: a one dimensional capillary fluid transport of a hardening adhesive 
O . through a single, straight vessel with diffusion of solvent through its walls and 

a mesoscopic scale based on topological characteristics of the vessel network, 
where results from the single vessel scale are mapped onto a periodic network. 
Given an initial amount of adhesive and applied bonding pressure, we calcu- 
late the portion of the filled structure. The model is applied to beech wood 
samples joined with three different types of adhesive (PUR, UF, PVAC) un- 

Sder various growth ring angles. We evaluate adhesive properties and bond line 
morphologies described in part I of this work. The model contains one free 
r^ , parameter that can be adjusted in order to fit the experimental data. 

Q ■ Keywords Adhesive • simulation • bond line • penetration model 

1 Introduction 

The most important principle in timber engineering to produce structural 
li*^ ■ wood components of constant quality, consists of cutting wood into smaller 

^s . pieces, selecting the best ones, and joining them again by adhesive bondings. 

• ■ What is known as a rather simple processing step becomes quite complicated, 

L^ . once we look in detail at the penetration of the hardening adhesive into the 

(^ ■ porous wood skeleton. Unfortunately the details of the adhesive penetration 
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can influence bond performance in multiple ways and the quality of the adhe- 
sive bonds determine the overall performance of structural parts [l]|2] . What 
compHcates studies of adhesive penetration is the interplay between pore space 
geometry and fluid transport, cell wall material and adhesive rheology and of 
course process parameters like amount of adhesive, growth ring orientation, 
and surface roughness, just to name a few [4]. While for soft wood predictions 
are rather simple, the micro-structure of hard woods complicates the prob- 
lem significantly, since adhesives can penetrate through the big vessel network 
deep into the wood structure jS] . In a previous work we explored the topolog- 
ical characteristics of the vessel network in beech (Fagus sylvatica L.) [6] and 
showed in part I of this work how the problem is dominated by flow through 
the vessel network. 

Adhesive penetration into hard wood was studied before [l3j|ll[71[8l[9] , al- 
though only experimentally or in descriptive form. For soft wood, the penetra- 
tion depth can be expressed by a simple trigonometric function, describing the 
filling of cut tracheids [10] . For hard wood however a model that characterizes 
the wood anatomy in order to predict the penetration depth and the amount 
of adhesive inside the structure is unknown. We construct an analitycal model 
based on the network properties and predict the adhesive penetration and the 
saturation of the vessel pore space. Our model has two-scales: the first scale 
describing the transport of a hardening adhesive through a single vessel in time 
due to an applied pressure and capillarity effects, and also with the possibility 
of constant diffusion of solvent through the vessel wall, what turns out to be 
important for some adhesives like PVAC or UF. When the viscosity increases 
by hardening and/or loss of solvent, the adhesive front slows down and finally 
stops. On the second, or network scale, the result for single vessels is embed- 
ded into a network model with identical topological properties like pore size 
distribution and connectivity that are characteristic for the vessel network of 
the respective wood. The model is compared with experiments where spec- 
imens are bonded with parallel longitudinal axes under varying growth ring 
angles using three different adhesive systems: PRF, UF, and PVAC. First we 
describe the rheological model of the adhesives, before we calculate the pene- 
tration into a single vessel with diffusion into the half space. Subsequently we 
discuss the network construction and the consideration of process parameters. 
With all model components at our hands, we finally compare the model with 
the experiments and discuss the results. 



2 Model Description 

Adhesive penetration is the result of an interplay of adhesive hardening, cap- 
illary penetration, and technological processing. In order to set up a model 
for adhesive penetration of hard wood, we have to combine several models 
in a hierarchical way. First we address bulk viscosity evolution of adhesives 
due to generic hardening mechanisms. On the fundamental level, we model 
the penetration of a fluid into one single, straight or wavy pipe. This model 



is enriched by diffusive transport of solvent through its wah. On the next hi- 
erarchic level, we project the fundamental model onto a network structure of 
perfectly aligned hard wood that represents the vessel network. Finally, we 
rotate the result of the vessel network penetration to consider the general sit- 
uation, where the adhesive surface is not necessarily aligned to the material 
orientation. We show how material parameters like porosity, hardening time 
or applied amount of adhesive will limit penetration. 



2.1 Modeling the hardening process 

The hardening process of various adhesives can be described by the tempo- 
ral evolution of the viscosity 77. Depending on the hardening type, different 
viscosity models need to be applied. For example reactive adhesives do not 
depend on the solvent concentration, while the viscosity evolution of solvent 
based adhesives strongly depends on solvent concentration C. In part I [11] 
of this work, we showed experimental viscosity measurements for UF, PVAC, 
and PUR. If solvent concentrations are important, like in the case of PVAC, 
the viscosity evolution can be expressed by 

7?(C,t) = r,g(C)[l + 7(C)exp(a(C)t)]exp(/?[l-C]) , (1) 

where rjg^ 7, a and /3 are parameters that depend on the adhesive type and 
the initial solvent concentration. For PVAC adhesive, we find 7=a=0, since 
the hardening process is mostly due to the loss of moisture and the initial 
viscosity only depends on the initial concentration. For PUR adhesive, the 
same expression can be used, however the concentration is kept constant during 
the process, expressed by /3=0 and constant 7^^, 7, a that only depend on the 
initial concentration. 

Unfortunately a whole class of adhesives, cannot be described by Eq. [H 
since their hardening process is more complex. For example the UF adhe- 
sive changes from liquid phase to gel phase during penetration, resulting in 
penetration arrest. The only active processes after this phase transition are 
the chemical curing reactions. Therefore the viscosity model should take into 
account the critical time when the phase transition occurs. Additionally, the 
concentration of the solvent changes in time due to the diffusion of the solvent 
into the cell wood structure. We propose the viscosity relation 

WO..) = .. -''^<^-7'-^'» , (2, 

Ci — t 

where g^i, 62, ai and ci are experimental parameters. Using the data from 
Ref. [11] we found (ii =6.985 x 10^ mPa-s^ and ai=10510s and variable param- 
eters (62 and ci) that depend on the initial solvent concentration. Note that ci 
describes the time when the penetration process finishes due to the liquid-gel 
transition. Using these two generic hardening models, we are able to describe 
the viscosity evolution of numerous adhesives. 




2.2 Single vessel penetration 

The fundamental scale is given by the capillary transport of a fluid charac- 
terized by its viscosity r^, inside a cylindrical pipe of radius R [12] with a 
penetration rate dl/dt that follows 

where ja = Pa -\- 2(j cos{0) / R with the applied pressure Pa^ the surface tension 
a and the contact angle between fluid and pipe wall 0. To obtain the penetrated 
distance l{t) we integrate 



(4) 



leading to a total fluid volume of Vp = 7TR'^l{t) inside the pipe. For reactive 
adhesives, whose hardening only depends on time, the integral can be found 
by combining Eqs. [T] and H] to 



R ^^ at-log(l + 7exp^^) V 
'^'^ = ^^[ a,,exp/3(i-C) ) • (^) 

This way we obtain the time dependent penetration distance in a straight sin- 
gle vessel, taking into account the applied pressure, the capillarity effects, and 
the reactive hardening process. Note that for adhesive types, whose viscosity 
changes when in contact with wood, Eq. H] can not be integrated so easily, 
since the viscosity depends also on the concentration that changes with time. 
Note that changes of the contact angle and furface tension of the adhesives 
with solvent concentration are not considered in this work. 

Since hard wood vessels are not straight, but weave tangentially around 
rays, the penetration distance needs to be modified. Here we simply describe 
vessels by the radius R^ wavelength A, and amplitude n (see Fig. [T]) of the 
oscillation in the z — ^/-plane in the parametrized form as, 

^2 ^ [y-ncos{kz)]^ ^ ^2 ^ .gx 

sec^ [arctan(n/c sin(/c2;))] 

where k = 27r/A. By integrating the vessel length ly along the z direction, we 
obtain the volume 



V{ly) = TVR^ly 



n'^P / sm{2kly) 



(7) 




unit cell 



Fig. 1 Vessel geometry and dimensions, vessel network, and unit cell. Shown are also the 
longitudinal (L), tangential (T), and radial (R) directions. 



2.3 Viscosity increase by diffusion 

Various adhesives contain solvents, whose concentration C in the mixture 
changes with time due to their diffusion into the cehular structure through 
vessel walls. To take this effect into account, we can write the solution of the 
diffusion equation in cylindrical z — r-coordinates as 



C(r, t) = — ;— exp 
^ ' ^ 2Dt ^ 



' ADt) 



(8) 



with the initial concentration of the solvent Co and its diffusivity across the 
cell wall D. The average diffusivity of the respective wood proved to be a good 
value. The mean value for the solvent concentration inside the vessel follows 



as 



C 



f 

Jo 



C{r,t)rdr = Co 



1 — exp 



' ADt 



(9) 



To obtain the complete equation for the evolution of the viscosity, we go back 
to Eq. [T] and insert the concentration evolution C into the respective concen- 
tration dependent parameters. Note that we do not consider the diffusion of 
low molecular parts of the adhesive. We also neglect the effect of swelling of 
the wood skeleton due to moisture changes, since the size of vessels is rather 
big compared to tracheids. 



2.4 Penetration into the network 



The adhesive penetration is dominated by the flow inside the vessel network, 
hence its topology determines the adhesive distribution. The network is formed 



by bundles of vessels that divide and weave around rays of various sizes. Inside 
the bundle, vessels interconnect by contact zones when touching each other 
and can also interchange positions J6l[T3]. Disorder in the network can only 
be considered through a numerical approach. In order to be able to derive 
an analytical model, we need to neglect disorder and use average topological 
network parameters. We build up a regular network using the average topo- 
logical parameters a and b for connectivity in tangential directions and c, d for 
the connectivity in radial direction. Fig. [T] shows the vessel network in three 
dimensions with the geometrical parameters a, 6, c and d. Note that a and b 
can be obtained from the size distribution of big and middle sized rays that 
are mainly responsible for the splitting and joining of the bundles of vessels. 
The parameters c and d however are more difficult to obtain. Basically the 
probability for radial network interconnections depends on the vessel density. 
We can therefore find a relation between the vessel density and the parameter 
c. d however will remain a free parameter for transport in radial direction. By 
separating the two geometric parameters a, b and c, d, we obtain anisotropic 
transport in the three principal directions, longitudinal L, radial i?, and tan- 
gential T (see Fig.[T]). Note that inclined samples with respect to the principal 
axis can be considered after rotation. 

To describe the penetration process of adhesives into wood, we have to 
define the bond line. The bondline is the whole region, where the adhesive 
can be found. This includes the pure adhesive between the two adherends 
and the area, where the adhesive has penetrated into the wood structure. The 
adherends are two pieces of wood which have been connected by the adhesive. 
In our case, we will focus on the zone where the adhesive layer and the adherent 
structure coexist. The procedure to obtain the maximum penetration depth is 
to calculate the penetration separately in each principal direction (tangential 
T, radial R^ and longitudinal L, shown in Fig.[T]) and then applying a rotation 
matrix to find the total penetration depth of the adhesive when the growth 
ring angle and the angle between vertical and longitudinal axis of the specimen 
are not zero. 

Since we employ a regular network, a unit cell can be used that consists of 
two single vessels with interconnections in the vertices and the center of the 
cell (see Fig. [1]). Consequently we can use Eq. Oof the single vessel to obtain 
the volume of the unit cell 



Vu = 2AbgT with gT = l^j^ (10) 



with the area A = ttR^. Because the unit cell can reproduce all the network, 
we can use it to simplify the calculation of the penetration depth on each 
direction and, using geometric properties, we can relate the adhesive volume 
inside the network, the network parameters and the penetration depth. 
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Fig. 2 Adhesive penetration into the vessel network with respective depths in tangential 
((5t), radial (^k), and longitudinal {bi^^ direction. Thick lines represent filled vessels. 



2.^.1 Penetration in the tangential direction 

To consider the penetration in the tangential direction (see Fig. [2j), we need 
to consider additionahy to the tangential waviness with wavelength b and 
amplitude a/4 the radial waviness with amplitude c/4 and wavelength d. Fig. [2] 
shows how the vessel network is filled by the adhesive. The path along one 
radial wave as function of the tangential coordinate xt is given by 



Sb = XT 



^idR 



1)1 



Attxt 






with gR = l 



C^TT^ 



8^2 • 

(11) 

The number of layers accessible from the bond line is defined by 2Lr/c with Lr 
the sample width. Summing over all accessible layers gives the total penetrated 
length as function of Xt 



2Lr 

Stotal = Sb[XT) 



(12) 



To obtain the total penetrated volume we have to calculate the number of 
unit cells Nul along Stotai and in longitudinal sample direction L^, this is 



NuL=^^^^, and then multiplied by the unit cell volume, Eg.lTUl 

Vt{xt)=N^lVu . (13) 



Assuming that the penetration of the adhesive is smaher than the total wave- 
length d, dsin{4jx/d) ^^^ Therefore using Eqs. M and [H Eq. M can be ex- 
pressed as 

Vt{xt) = XT — LnLLgrgn • (14) 

ac 

When the adhesive stops to penetrate, this volume becomes the maximum 
volumen V inside the structure and the tangential coordinate xt transforms 
in the maximum penetration depth St, 

AALnLLgrgn 

2.4-2 Penetration in the radial direction 

Following the idea of calculating the adhesive penetration in each principal 
direction, the next step is to obtain the penetration depth 5r when the adhesive 
penetrates only in radial direction. Fig. [2] illustrates the penetration of the 
adhesive in order to relate the volume with the network parameter and the 
penetration depth. We analogously count the volume occupied by vessels as 
function of the radial coordinate xr. Again we calculate the total length Stotai 
of the radial wave but now as function oixR by Stotai^'^XR/ csi){Lt) and obtain 
the number of unit cells 

NuR[xr) = — . (16) 

cab 

The total volume occupied is given by multiplying the number of unit cells 
NuR (Eq. [HI) by the volume K from Eq. [JOl 

4A 

Vt{xr) = VuNuR = XR — LLSb{LT)gT • (17) 

ac 

As before, we must now compare Vt{xr) with the volume occupied by the 
adhesive V with the maximum penetration depth in the radial direction xr = 
Sr to obtain 

4A Vac 

V = 5r — LLSb{LT)gT , and 5r = — — -— - — . (18) 

ac 4ALLSb{LT)gT 

Finally, we can insert the penetration path from Eq. [TT] and obtain 

Sn = . '^, . (19) 



AALlLt 



l^^{gR-l) (l 



AttLt 



gr 



2.4-3 Penetration in the longitudinal direction 

We consider now the penetration only in longitudinal direction. Fig. [2] shows 
that the adhesive penetration Sl is basicahy along the vessels. This value is 
found again by calculating the number of total unit cells, but now in the plane 
LtLr, namely 

iV„. = ^^^(M . (20) 

c a 

Multiplying Nul with the occupied volume V{xl) of the adhesive for each 
vessel as function of xl (Eq. [7]), and taking into account that the penetration 

Vt{xl) = xl—Lrs,{Lt) {2gT - 1) • (21) 

ac 

Again, comparing this volume with the adhesive volume, using Eq. [TTJ we can 
write the maximum penetration depth as 



Vac 

ol = 



4ALrLt 



i + H..-i)(i + ^^4g^ 



(22) 



{2gT - 1) 



Finally we obtained the maximum penetration depth in the three principal 
directions fEqs. ll5|19|22p . We can introduce the porosity e of the wood which 
can be extracted easily from experimental data [6 . Expressing the penetration 
depth in terms of porosity also simplifies the model verification. The number 
of vessels Ny in the plane LtLr equals Ny = 2- Nul- The porosity is therefore 



NyA 4 A 



LrLt ac 



-^(9R-l)\l + 



(23) 



Since porosity is a mean value, we can neglect the periodic part on the right 
hand of the Eq. [23l Inserting e into Eqs. [151 [HI and [22l the maximum pene- 
tration depths become 

^^^ V ^^^^ V ^^^ ^ V{9R^l) ^ (24) 

elLLrgr ' ^LrLt {2gT - 1) ' 2eLRLLgTgR 

2.4-4 Limitation due to the total amount of applied adhesive 

Up to now, we calculated the penetration of an infinite amount of non-hardening 
fluid. However the amount of applied adhesive and the penetration time due 
to hardening are both limited. Therefore the volume V needs to be calcu- 
lated considering these limitations. Both limitations will lead to different pen- 
etrated volumes but only the smaller one has a physical meaning. To calculate 
the volume with penetration of hardening adhesives, we need to treat the 
L, i?, T— directions separately. 



10 

To consider tangential penetration we employ Eq. [2l] and apply adhesive 
only on the RL plane. From there, the adhesive can penetrate two channels 
with radius R per unit cell, and considering the number of unit cells on this 
face NuT^ the volume penetrated after the hardening process, using Eq. |H is 

Vrl = Vp^^ . (25) 

Inserting Eq. [25] into Eq.lMland using Eq.[23l we obtain the penetration depth 
with hardening Srh as 

^Th = T • (26) 

grgn 

For the radial penetration^ the penetrated volume is given by. 

By inserting Vlt into Eq. [24] and taking the mean value of the periodic term, 
we obtain 

^^'^ = "b {gn + DsT ■ ^''^ 

Finally, for the longitudinal penetration^ J^h, following a similar procedure the 
penetration depth becomes 

m 



JLh 



{2gT - 1) 



(29) 



These values determine the maximum penetration depth that the adhesive 
can reach until becoming solid. However it is possible, that not enough adhesive 
is available, and penetration stops before. Using the available adhesive volume 
V in Eqs. [2H the penetration depths Sr, St, and Sl can be calculated and 
compared to the hardening ones {Sr^, Srh, ^Lh), e.g. if SrKSr^, to obtain the 
limiting case. 

2.5 Penetration depth for an arbitrary orientation 

In order to apply our model to real situations, we must have a way to consider 
an orientation of the adhesive application surface that deviates from the wood 
material system. Therefore we need to calculate the global penetration depth 
Sy and Sh as function of 6r, 6l, ^t, and Sri^^ ^^h, ^Th, respectively. We can 
apply a rotation matrix with the growth ring angle ip and the angle between 
the vertical axis and the longitudinal axis of specimen. Assuming that the 
adhesive is always applied on the y — z plane, we apply two rotations in the 
principal coordinate system, one in the radial direction and the other in the 
longitudinal direction via the rotation matrix 

' sin(?/;) — cos{6) cos('^) — cos('^) sin(6>) ^ 
M= I cos(V^) cos(l9) sin(V;) sin(l9) sin(V^) | . (30) 

- sin((9) cos(i9) 
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Note that Eqs. [2l] give a dependence of the penetration depths on the 
apphcation areas Lr - Lt, Ll ■ Lt, Ll ■ Lr. We define a penetration vector 
6 = Sy^, where S_ is oriented normal to the adhesive surface. S is given by 

V f 1 lign^l) 1 



S=-l-^A25^^^ \ (31) 

^ \9t grgn 2gT-lJ 

in the principal coordinate system {R — T — L). Applying the rotation matr 
ix M to the vector 6_, the x component gives the maximum penetration depth 

At Al Ar 

Sv = —. — cos(^) sin('0) + — — sin(^) sin('0) + — — cos('0), (32) 

^ad ^ad ^ad 

with the area Aad of the surface where the adhesive is applied. We can directly 
apply the rotation matrix to the vector 6h={SRh^ ^Th, ^Lh) and find, 

^h = ^Th cos(6>) sin(?/^) + 5Lh sin(6>) sin(V^) + Sri, cos(?/^). (33) 

With these derivations, we complete the geometric and dynamical descrip- 
tion of our model. The information about the dynamics of the adhesive is 
included in the length l{t) according to Eq. HI Finally, our maximum penetra- 
tion depth with solvent diffusion can be calculated using Eq. [33l by replacing 
the concentration function in Eq. [T]for the respective adhesives. In a next step 
we will apply the model to experiments described in the first part of this paper 



3 Application of the Model 

Using synchrotron radiation X-ray tomographic microscopy (SRXTM) and 
digital image analysis, we extracted bond lines from beech wood samples, that 
were bonded with PUR, UF, and PVAC adhesives of different viscosity under 
growth ring angles ranging from 0° to 90° in 15° steps [llj. Since our model is 
periodic, we will calculate the maximum penetration depths for various situa- 
tions. The procedure is as follows: First we calculate the penetration distance 
l{t) of adhesive inside a single vessel. Note that for PUR the calculation is with- 
out time dependence of the concentration using Eq. O while for PVAC and 
UF adhesives with time dependence additionally Eq. [9]is used. The porosity e 
and mean radius of the vessel R are taken from an earlier SRXTM study [6 as 
i?=28.03/im and porosity e=0.34. For all samples the mean applied pressure 
was PA=0.7MPa[ll . Literature values of the surface tension a |lH|T5[|16], for 
the three types of adhesives, are not large enough to compete with the applied 
pressure term in Eqs. [3l leading to negligible capillarity effects, this means 
Ii=Pa' The parameters for the viscosity rjg^ a, j3 and 7 are taken from part I 

mj. 

— For PUR adhesive, we choose the concentration C — 0.71, and the param- 
eters r]g = 4911mPa-s, 7=9. 74x 10-^ a=0.0028s-\ f3=0. The diffusion of 
solvent is not relevant. With these quantities the length for PUR adhesive 
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Fig. 3 Dependence of the viscosity parameters 62 and ci with the solvent concentration 
for UF adhesive. The dots denote the experimental data and the solid line the exponential 
fit. Experimental values for UF are ui = 4.542, z/2 = 0.5379,^/3 = 358. 8,1/4 = 3.482. 



K^)puR is calculated using Eq.[5]to /(t)puR=0.304m. This value seems huge 
at first sight, however it related to the path along the waving vessels that 
can easily reach lengths of 0.5m and above. 

— In the case of PVAC adhesive, the parameters are Co =0.49, 7^^=0.001859 
mPa-s, 7=0, Qf=Os~^ and /3=29.64. The diffusivity of the solvent (water) 
for the samples is taken as D=3.0xl0~^ m^ s~^ [HUE,, and using Eq. |H 
we find a significantly lower value /(t)pvAC=0.7mm. 

— For UF we include the viscosity parameters 62 and ci of Eq. [2] that change 
with the solvent concentration. We use the experimental viscosity data 
from Ref.[TT] and fit it with analytical curves (see Fig.[3|) to determine the 
concentration dependence of the viscosity parameters 62 and ci . After the 
identification of the right values for 62 and ci , we can integrate Eq. |4] and 
obtain a vessel penetration depth of /(t)uF=l-lnim. 

The parameters a and b can be determined experimentally using image 
processing (see Ref.[TT]). In our case we measured the area and the eccentricity 
of segmented rays and averaged over several samples, and obtained values for 
a=0. 156mm and 6= 1.5 74mm. To eliminate variations due to the year ring 
structure, we used an average porosity of e=0.34. The cylindric sample size had 
10mm height and 3mm diameter, leading to an adhesive area of Aad=^Om.m^. 
As described in Ref.[Tl , the quantity of applied adhesive was around 200g/m^ 
for all adhesives. 

We compare the maximum penetration depth for samples with different 
growth ring ip and grain angles (see Figs. l4l5p . We fit the parameter d to ob- 
tain ^^ = 1 , what can be interpreted as a lower probability of interconnection 
in radial than in tangential direction. 

— For PUR adhesive, we choose a sample with angles, ?/^=50° and 6^=3.85°. 
If we calculate the maximum penetration depth using Eq. [32] we obtain 
a value of (5/,, =38. 5mm with hardening as limiting factor, however using 
the volume limitation with Eq. [33] we obtain (5y=848/im. Therefore we 
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Fig. 4 Bond lines with PUR adhesive in beech wood samples. The maximum penetration 
depth predicted by the model for all samples with growth ring angles -0 and grain angle 
is shown by the white lines. 



can conclude that all the adhesive penetrated before hardening took place, 
leaving a starved bond line behind (see Fig. |4]). Note that the adhesive 
penetrates both wood pieces, but significantly deeper into on the applica- 
tion side (right side of samples in Fig. S]). Therefore we have to take an 
average value of (5=424/i?7i showing good agreement with the experimental 
data. To test the model for other orientations ip^O^we choose a sample with 
ip=l.l° and ^=3.6°. This means we use the previous calculation but apply 
a new rotation matrix M. We obtain a penetration depth of Jy=596/im, 
^=298/im and (5/i=1.32mm. In Fig. |3]the quality of the analytical predic- 
tion is shown. We repeat this for angles '0=27.9° and 6>=3.8° and obtain 
the penetration depths 5v=Slljiim^ J=405.5/im and (5/^=23. 6mm (compare 
Fig. m. These tests show that our model is a good approximation for the 
beech wood structure and therefore we fix the network parameters a, 6, c, d 
for further calculations. 

— We exemplify the penetration of PVAC using a sample oriented at angles 
?/^=46.1° and ^=3.4° and calculate the maximum penetration depths from 
Eqs.[32land[33l We find (5y=851/im, leading to 5=425.5 jam and J^=80/im. 
Therefore the maximum penetration depth is limited by the hardening 
process. In Fig. [5l we show that almost all the adhesive remains in the 
bond line with only a small quantity of adhesive inside the vessel network. 

— For UF we repeat the same procedure as before on a sample with orienta- 
tion angles '0=37.2° and ^=4.3°. We find that the penetration depths are 
^y=856/im, J=428/im, and J/^=141/im and again the penetration of the 
adhesive is limited by adhesive hardening. Fig. [5] shows the sample with 
the predicted penetration depth, exhibiting excellent agreements between 
the analytical prediction and the experiments. 

To study the dependence of maximum penetration depth on the growth ring 
angle, we take the values for UF and keep all parameters fixed, except the 
growth ring angle. In Fig. [6] we see the two limiting conditions for the pene- 
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Fig. 5 Bond lines of PVAC and UF adhesive in beech wood with maximum predicted 
penetration depth for samples with the orientation angle ip and 0. All dimensions are given 
in /im. 



tration depth. The penetration depth is an increasing function of the growth 
ring angle for the hardening hmitation case, and we observe a distinct maxi- 
mum at approximately 48° in the case when the maximum available volume 
is the limitation. This observation is in agreement with PUR adhesive that 
fulfills the volume limiting condition, as demonstrated in part I of this work 
pT] . This result shows that even though we reduce the wood anatomy to a 
homogeneous, regular network, adhesive transport, the beech wood seems well 
described by the model and for a desired penetration depth, the model can 
predict the optimal growth ring angle of the samples. 
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Fig. 6 Dependence of the penetration depth on the growth ring angle. 



Our model can also be used to design new adhesives with optimized prop- 
erties, like reactivity, if an ideal penetration depth is to be reached. Fig. [71 
shows the maximum penetration depth for a wide range of adhesive parame- 
ters i^i, ..., i^4 from Fig.[3l We show this in four plots combining two parameters. 
Horizontal planes represent the case where all available adhesive is inside the 
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vessel structure, while the curved surfaces show the penetration limit due to 
adhesive hardening. The intersection line (see Fig. [7j) separates regions with 
complete penetration from those, where penetration is limited by adhesive 
hardening. Therefore, Fig. [71 allows to choose a pair of reactivity parameters 
in order to obtain a desired penetration depth. The model can also be used to 
minimize solvent concentration and amount of applied adhesive for a required 
penetration depth. Fig. [8] illustrates the maximum penetration depth as func- 
tion of the solvent concentration and the total amount of applied adhesive. 
The solid lines represent the proportions between solvent concentration and 
the total applied volume of adhesive which give the same penetration depth. 




0.0 
1000 



Fig. 7 Penetration depth 6 in mm as function of adhesive parameters from Fig. [3] for UF 
and growth ring angle of 45°. The experimentally obtained value for UF is marked with the 
white dot. (Color version online) 



4 Discussions and Conclusions 



We presented an analytical model for the prediction of the penetration depth 
of adhesives, paint or hardening fluids in general, into the beech wood struc- 
ture. Since we focused on hard wood, the pore space of wood is formed by 
a network of interconnected vessels. The network is characterized by param- 
eters that are related to amplitude and wavelength of the oscillating vessels 
in tangential and radial direction and the porosity of wood originating from 
the vessel network. We compared the model to experiments and found good 
agreement for various adhesives. Therefore, even though we reduce the wood 
anatomy to a homogeneous, regular network, adhesive transport for the much 
more disordered, complex pore space of real beech wood seems well described. 
The analytical model considers generic types of adhesive hardening. How- 
ever if other special fluids are to be considered, only the viscosity dependence 
on concentration and time need to be known. We applied the model to three 
major types of adhesive which are PUR, UF, and PVAC and compared the 
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Fig. 8 Penetration depth 6 in mm as function of solvent concentration and total amount 
of applied adhesive for UF and growth ring angle of 45°. The solid lines represent the fixed 6 
and the dashed line divides the surface into two regions, the right one when the penetration 
is limited by the hardening process and the left one, when the total volume of applied 
adhesive is the limiting factor. (Color version online) 



respective penetration depths. For adhesives whose hardening process depends 
on the change of concentration, we include a description for solvent concentra- 
tion diffusion inside the beech wood. This approach makes the model appli- 
cable for adhesives like UF and PVAC. Penetration is limited by two things: 
the penetration due to the applied pressure that is arrested by hardening pro- 
cesses, and the total amount of applied adhesive that is available to penetrate 
into the vessel network. The smaller penetration depth is the limiting one. By 
comparing the model with the experimental data we showed that it is possible 
to model the maximum penetration depth for three different used adhesives, 
namely PUR, PVAC, and UF. 

Our model is sufficiently simple to allow for a broad applicability. By de- 
termining the morphological and rheological parameters, it can be applied to 
a wide range of wood species and to fluids with various hardening kinematics 
to predict the penetration depth of these fluids into porous structures, when 
transport is dominated by capillarity. 
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